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Screenings of the quark-gluon plasma in electric and magnetic sectors are studied on the basis 
of generalized Polyakov-line correlation functions in lattice QCD simulations with two flavors of 
improved Wilson quarks. Using the Euclidean-time reflection {TZ) and the charge conjugation (C), 
electric and magnetic screening masses are extracted in a gauge invariant manner. Long distance 
behavior of the standard Polyakov-line correlation in the quark-gluon plasma is found to be dictated 
by the magnetic screening. Also, ratio of the two screening masses agrees with that obtained from 
the dimensionally-reduced effective fleld theory and the A/" = 4 supersymmetric Yang- Mills theory. 
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I. INTRODUCTION 

Screenings inside the quark-gluon plasma (QGP) at 
high temperature are dictated by thermal fluctuations of 
quarks and gluons, and are characterized by the electric 
and magnetic screening masses [l| . From the phenomeno- 
logical point of view, screening of color charge affects the 
force between heavy-quarks and determines the fate of 
heavy-quark bound states such as J/ ip and T inside hot 
QCD matter produced in relativistic heavy-ion collisions 
at RHIC and LHC [111. 

Non-perturbative determinations of the screening 
masses have been attempted from various approaches in 
the past: Use of Polyakov-line correlations with appropri- 
ate projection in color space [343|) use of correlations of 
spatially local operators classified in terms of Euclidean 
time-reflection iSj], use of dimensionally-reduced effective 
field theory [l^ and a direct lattice evaluation of the 
gluon propagator under gauge fixing fTl|. The screening 
masses of the A/" = 4 supersymmetric Yang-Mills theory 
for large 't Hooft coupling have been also analysed on the 
basis of the AdS/CFT correspondence [HI- 

In this paper, we extract electric and magnetic screen- 
ing masses from gauge-invariant Polyakov-line correla- 
tions in two-flavor lattice QCD simulations (a prelimi- 
nary account has been reported in [11). According to the 
symmetry properties under the Euclidean time-reflection 
TZ and the charge conjugation C [1], we decompose the 
standard Polyakov-line operator Q, into four independent 
types, r^Mit and risi- Here M and E stand for the 7?.-even 
magnetic sector and the TZ-odd electric sector, respec- 
tively, while ± stands for the even or odd under C. Then 
the magnetic (electric) screening mass is extracted from 
the correlation of Trr2]vi+ (TrilE- ) ■ In our simulations on 
a 16'^ X 4 lattice, we employ two-flavor lattice QCD with 
improved Wilson quarks coupled to an RG-improved glue 



(Iwasaki gauge action) for the temperatures T/Tpc ~ 1- 
4 with Tpc being the pseudocritical temperature. We 
take the quark masses corresponding to nip^/m^ =0.65 
and 0.80 where rnps (m-v) is pseudoscalar (vector) meson 
mass at T = 0. 

This paper is organized as follows: In Sec. [Hi we dis- 
cuss the properties of the Polyakov-line operators under 
TZ and C. Numerical simulations and the results of the 
screening masses are shown in Sec. IIIII Comparison of 
our results with those obtained from the dimensionally- 
reduced effective field theory and the TV = 4 supersym- 
metric Yang-Mills theory are also given. Section IIVI is 
devoted to summary and concluding remarks. 



II. GENERALIZED POLYAKOV-LINE 
CORRELATIONS 

We start with the standard Polyakov-line operator 

f2(x) = Pexp ig 



f / dTAi{T,:x.) 
Jo 



(1) 



Its gauge invariant correlation function reads 
Cn{r,T) = (Trl7t(x)Trl](y)} - \{Tin)\^, 



(2) 



where r = |x — y| and we define Trfi = X]a=i ' 
with a being color indices. The screening mass is de- 
fined from the exponential fall-off of this gauge invariant 
correlation, 



Cnir,T) j,,iT)- 



rT 



(3) 



Since the magnetic gluon is less screened at finite T, the 
dominant contribution to m^^ (T) would come from non- 
perturbative magnetic sector of QCD as pointed out in 
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In order to extract the electric and magnetic screening 
masses separately, we classify the Polyakov-line operator 
into different classes on the basis of the symmetries under 
TZ and C which are good symmetries of QCD at zero 
chemical potential |8j]. Under TZ and C, the gluon fields 
transform as, 

yl,(r,x) ^ A,(-T,x), A4(t,x) ^ -A4(-r,x),(4) 
A^(t,x) 4 -A;(r,x). (5) 

We call an operator magnetic (electric) if it is even (odd) 
under TZ. Under these transformations, the standard 
Polyakov-line operator transforms as 

n^n*. (6) 

Then, we can define magnetic (electric) operator which 
is 7?.-even (7^-odd) as 

nM = l{n + n^), nE = ^{n-n^). (?) 

Furthermore, they can be decomposed into C-even and 
C-odd operators as 

^M±^^{nM±n*M) r!E± - i(f^E±r!^). (8) 

Using these operators and noting the fact that 
Trf^M- = Trr2E+ — 0, we can define two generalized 
gauge-invariant Polyakov-line correlation functions, 

CM+{r,T) = (Trr!M+(x)Trr!M+(y))-K™)|2, (9) 
CE-{r,T) = (Trr!E-(x)Trl]E^(y)). (10) 

Note that Tri7M+ (Tri^E-) is nothing but the real (imag- 
inary) part of Tiil, and that (Trfi) is real due to the C 
symmetry. 

The electric and magnetic screening masses can be de- 
fined from the above gauge invariant correlation func- 
tions through their long distance behavior: 

g-mM+(T)r 

CM+{r,T) > 7,,^(T) , (11) 

~mE-(T)r 

CE-(r,r) > 7,_(T) , (12) 

r— >-oo Tl 

Notice that the standard Polyakov-line correlation 
function and the above generalized correlation functions 
are simply related as 

Cn{r,T) = CM+{r,T)-CE-{r,T). (13) 

Thus the separate determination of Cm+ (r, T) and 
CE-{r,T) on the lattice enables us to study the rela- 
tive importance of the magnetic and electric sector in a 
non-perturbative manner. 



III. LATTICE SIMULATIONS 
A. Lattice setup 

We employ a renormalization group improved gauge 
action Sg and a clover improved Wilson quark action 
with two fiavors Sq: 

4 4 

/=1,2 x,v 

where /3 = 6/g^, ci = -0.331, cq = 1 - Scj, and WJ^^"" 
is a n X TO rectangular shaped Wilson loop and 

+ (1 + lt^Ul^^6x^y+[i} - SxyCswK E <^f^iyFtJ.u- (16) 

Here K is the hopping parameter and F^^, is the lattice 
field strength, F^^ = l/8i{f^^ - f^J, with /^^ being 
the standard clover-shaped combination of gauge links. 
For the clover coefficient csw^ we adopt the mean field 
value using VF^^^ which was calculated in the one- loop 
perturbation theory, 

csw = (W^^^i)"3/4 ^ (1 _ 0.8412/3-i)-3/4. (17) 

Our simulations are performed on a lattice with a size of 
X Nt — 16'^ X 4 along lines of constant physics, i.e. 
lines of constant TOpg /m^ (the ratio of pseudoscalar and 
vector meson masses) at T = in the space of simulation 
parameters. Details on the lines of constant physics and 
the phase diagram for iVj = 4 with the same actions as 
above are given in Refs. 0, [l3|- We take two values, 
TOpg/wy — 0.65 and 0.80, with the temperature range 
of T/Tpc - 1.0-4.0 (10 points) and 1.0-3.0 (7 points), 
respectively, where Tpc is the pseudocritical temperature 
along the line of constant physics 0, [3- The number 
of trajectories for each run after thermalization is 5000- 
6000, and we measure physical quantities at every 10 
trajectories. The statistical errors are estimated by a 
jackknife method with a bin size of 100 trajectories. 

B. Screening masses 

Results of the generalized Polyakov-line correlations 
Cm+(e-)(''j 2^) are shown in Fig. [1] for m^^/m^ — 0.65 
(upper panel) and 0.80 (lower panel). These figures show 
that (i) the magnetic correlation and the electric corre- 
lation have an opposite sign, and (ii) the magnetic cor- 
relation has larger magnitude and longer range than the 
electric correlation at long distances. The latter implies 
that the standard Polyakov-line correlation is dominated 
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FIG. 1: Results of CM+{r,T) and CE-{r,T) for several tem- 
peratures as a function of rTpc at ra-p^/m^ — 0.65 (upper 
panel) and 0.80 (lower panel). 



FIG. 2: Results of CM+(r,T) (upper panel) and ~CE-{r,T) 
(lower panel) in a logarithmic scale as a function of rT at 
mpg/my = 0.80. 



by the magnetic sector. In Fig. [2] we show CM+{r,T) 
(upper panel) and — CE-(r, T) (lower panel) in a loga- 
rithmic scale as a function of rT for m-p^/m^ = 0.80. 
These correlators behave linearly and are scaled well by 
rT at high temperature, so that the screening masses can 
be extracted by fitting the correlators using the Yukawa 
form defined in Eqs. ((TTl) and ([T2|) . 

We fit these correlators by Eqs. pT|) and in 
the same interval of 0.5 < rT < 1.0 by minimizing 
X^/^dof-^ Upper (lower) panel of Fig. [3] shows results 
of screening masses in the magnetic and electric sector 
at mpg/m^ = 0.65 (0.80) as a function of temperature. 
Numerical results are summarized in Tab. U and |TT] for 
TTipg/my = 0.65 and 0.80, respectively. The screening 
mass obtained from the standard Polyakov-line correla- 
tion function m^{T) using Eq.([3]) is also shown. As ex- 
pected, the magnetic sector has longer range than the 
electric sector {mM+{T) is smaller than mE-(T)) and 



We study the fit range dependence of the results, and find that 
the magnitude of systematic error due to the fit range is smaller 
than or comparable to the statistical errors at T ^ 1.2Tpc. 



also the standard screening mass is dominated by the 
magnetic mass {m^(T) ~ mu+iT)). 

C. Comparison to other approaches 

Let us compare our results with the predictions by the 
dimensionally-reduced effective field theory (3D-EFT) 
and the Af = 4 supersymmetric Yang-Mills theory 
(SYM). 

In the 3D-EFT approach [I], the screening masses m 
various channels with the quantum numbers J^*^ are 
calculated: Here P is a parity in a two-dimensional 
plane perpendicular to x-y, and J is an angular mo- 
mentum in the two-dimensional plane. Since the lowest 
masses for {7l,C) = (+,+) and {TZ,C) — {—,—) are in 
(J, V) ~ (0, +) channels, we extract ttim-i- from the O^"*" 
channel and rriE- from the 0^~ channel. Tables 6 and 
8 of Ref. P for the lattice coupling 13 = 21, lattice size 
L = 30, iV/ = 2 and T ~ 2Aj^ lead to 

3D-EFT {Nf = 2) : (18) 
m„jr = 3.96(5), TO,_/r = 7.01(10), 
'^E-/™M+ - 1-77. 
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FIG. 3: Results of magnetic and electric screening mass 
m-^j_/T and m-^_/T together with the screening mass ob- 
tained from the standard Polyakov-line correlation m^/T as 
a function of temperature at m-p^/m^ — 0.65 (upper panel) 
and 0.80 (lower panel). 



TABLE L Results of screening masses at ra-pg/m^ — 0.65. 
First parentheses show statistical errors, while the second 
parentheses show systematic errors calculated from difference 
between screening masses in the range of 0.5 < rT < 1.0 and 
V6/4 <rT < 1.0. 
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II: Results of screening masses at m-pg/m^ 


errors calculated by the same procedure as 
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The screening masses of the Polyakov-line correlation 
functions in TV = 4 supersymmetric Yang-Mills theory in 
the limit of large Nc and large 't Hooft coupling (A = 
g^Nc) were calculated by AdS/CFT correspondence [l^ . 
Under the same identification in (J, V) — (0, -f ) channel 
as discussed above, Tab. 1 of Ref. [I4I leads to 

A/" = 4 SYM : (19) 

m„jr = 7.34, m^_/r = 16.05, 

"^E-/™M+ = 2.19. 

Results of our screening masses and their ratio for 
1.67 < T/Tpc < 3.22 with mpg/m^ = 0.65 in Tab. U 
indicate that 

4D-lattice QCD {Nf = 2) : (20) 
TO„JT- 5.8(2), TO,_/T = 13.0(11), 
"iE-/"tM+ = 2.3(3). 

In all three cases above, we have an inequality < 
TOj,_ so that the magnetic sector dominates at large dis- 
tances. Since we cannot compare the absolute magni- 
tude of the screening masses in QCD and that in A/" = 4 



SYM due to different number of degrees of freedom, we 
take a ratio, ra-^_ /m^j^ (Table Hill) . and make a compar- 
ison for the three cases in Fig. We find that the ratio 
agrees well with each other for the temperature range of 
1.5 < T/Tpc < 3. 



IV. SUMMARY 

We investigated the screening properties of the quark- 
gluon plasma from the Polyakov-line correlation func- 
tions classified under the Euclidean-time refiection {TZ) 
symmetry and the charge conjugation (C) symmetry. The 
Polyakov-line correlators with the quantum numbers of 
(TZjC) = (even, and (odd, — ) are defined in gauge in- 
variant form. 

From the lattice simulations of Nf — 2 QCD above 
Tpc with the RG-improved gluon action and the clover- 
improved Wilson quark action, we extracted the mag- 
netic screening mass m^,^ and the electric screening mass 
TOj,_ , and found that m-^_ /m^^ = 2.1(2) for 1 < Tpc < 3. 
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TABLE III: Results of the screening ratio /m^,^ in Nf — 
2 lattice QCD simulations. 
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(3D-EFT) 01 and Af = 4 supersymmetric Yang-Mills theory 
(SYM) 



Therefore, the standard Polyakov-line correlation func- 
tions at long distance is dictated by the magnetic screen- 
ing. Our electric-magnetic ratio is consistent with those 
obtained from the dimensionally-reduced effective field 
theory (3D-EFT) and the A/" = 4 supersymmetric Yang- 
Mills theory (SYM). 

In the present work, wc have not attempted to make 
projections of the Polyakov-line operator 57 to the opera- 
tors with definite J (two-dimensional angular momen- 
tum) and V (two-dimensional parity). It would be a 
future task to make such projections, so that we can 
compare the 4D-lattice resuhs with 3D-EFT and SYM 
in more details. In Ref. [13], we extracted four screen- 
ing masses (X = M-I-, M— , E+, E— ) from gauge-fixed 
correlation functions Cx{r,T) = (Tr[rix(x)rix(y)]}; we 
found that the long distance behavior of the color-singlet 
Polyakov-line correlation functions is dictated by m^^ 
which is smaller than and 77ij,_ for T > l.STpc. Fur- 
ther studies are needed to have physical interpretation of 
this result. 
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